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GEOMAGNETIC ACTIVITY: 

DEPENDENCE ON SOLAR WIND PARAMETERS 

by Leif Svalgaard 

Institute for Plasma Research 
Stanford University 

1. THE MAGNETOSPHERE 

The magnetized collisionless solar wind plasma con- 
fines the magnetic field of the earth (and other solar 
system bodies) to a region around the planet called a 
magnetosphere. Alternatively we could say that the 
planetary magnetic field excludes the solar wind from 
the planetary environment. The confinement of the 
field or the exclusion of the solar wind plasma is not 
perfect, however. Due to particle gyration, the two 
plasma regimes overlap slightly thereby allowing signi- 
ficant mutual interaction. In spite of extensive ef- 
forts there is still no satisfactory physical theory 
of the interaction between the solar wind and' the 
magnetosphere. In the broadest terms it seems that 
the magnetosphere may be described as a resistive - 
and therefore dissipative - element in an electrical 
circuit drawing its power from the kinetic energy of 
the streaming solar wind plasma. The tendency of a 
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plasma to react in highly dissipative very often 
explosive - ways to any changes in boundary conditions 
is well-known from laboratory experiments and is also 
amply manifested in the magnetosphere as it continu- 
ously readjusts itself to the everchanging solar wind. 

In a sense the magnetosphere is a continuing 
plasma physics experiment. The plasma parameters as 
well as the basic geometry of the experiment vary on 
many different time scales. Monitoring of the experi- 
ment is provided *by observations of geomagnetic and 
auroral activity, and recently also by measuring di- 
rectly the relevant external plasma parameters. 

It has been established that geomagnetic activity 
indeed can be used as a measure of solar wind pa- 
rameters, and we are now in the position of being 
able to "calibrate" the experiment. It has long 
been recognized that "certain restricted areas of 
the Sun's surface are responsible for terrestrial 
magnetic disturbances" (Bartels, 1932) which then 
"yield supplementary independent information about 
solar conditions". 

Penetration of a streaming plasma into a region of 
strong magnetic field may be considered to take place 
in a way lying between the following two extremes: 

1) The plasma is completely diamagnetic and excludes 
the field from its interior by flowing around the 
field region, or 2) The plasma remains non-diamag- 
netic as it encounters the magnetic field and crosses 
it by means. of an electric polarization and corre- 
sponding E*B drift. The first viewpoint has been 
rather successfully used to account for the- basic ex- 
istence as well as for the approximate size and shape 
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of the magnetosphere, but fails in providing a 
physical basis for geomagnetic activity and disturb- 
ances. The second extreme seems to account for the, 
behavior of laboratory plasma streams (e.g., Baker 
and Hammel, 1965) j and may be applicable to the mag- 
netosphere as well (Eastman et at, 3 1976). 

Consider a plasma stream moving. into a strong uni- 
form magnetic field normal to the plasma velocity 
vector. A schematic cross-section is .shown in Figure 
1. The effect of the Lorentz V*B force is canceled 
by the action of an electric field E = - V*B that, if 
it is not already present, will be produced by a slight 
charge separation in the plasma caused by the Lorentz 
force deflecting ions and electrons in opposite di- 
rections. Now the plasma can move across the magnet- 
ic field due to the well-known ExB drift. Electro- 
static repulsion will tend to spread the polarization 
charge layers laterally along the magnetic field lines. 
If the field lines pass through a good conductor as 
shown in Figure 1, the polarization charge can be neu- 
tralized and the plasma stream will consequently be 
stopped or retarded. The depolarizing current along 
the field lines can thus serve to transfer momentum 
from the plasma stream to the conductor. 

Figure 2 shows this model applied to the magneto- 
sphere. As the solar wind plasma streams over the 
polar regions the dawn side of the magnetosphere be- 
comes positively charged and the dusk side becomes 
negatively charged by the VxB force. Depolarizing 
currents along high-latitude magnetic field lines 
close through the polar ionosphere hence imparting 
some of the solar wind momentum to the ionospheric 
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CONDUCTOR 


Figure 1. Schematic of a plasma stream 
moving into a transverse magnetic field B. 

The polarization electric field E = - V*B 
follows from a charge separation in the 
plasma. The charge layers can spread along 
magnetic field lines (the current i) and 
may be neutralized if a good conductor con- 
nects the two spreading charge layers. 

and magnetospher ic plasma and producing a magneto- 
spheric boundary layer of retarded solar wind plasma. 
The geomagnetic field lines through the boundary 
layer become stretched downstream forming the magneto- 
tail. Thus, the magnetosphere acts like a magneto- 
hydrodynamic generator, converting kinetic energy of 
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moving plasma into electrical energy which then ul- 
timately is dissipated within the magnetosphere or 
in the boundary layer. When the northern and southern 
parts of the boundary layer meet downstream of the 
earth, plasma regions with oppositely directed magnet- 
ic field lines are brought together and may reconnect 
explosively. The resulting "magnetospheric fireball" 
(Frank et al. 3 1976) is a source of energetic par- 
ticles which, if precipitated into the earth's atmo- 
sphere, produce the often spectacular auroral dis- 
plays associated with geomagnetic disturbances. 



Figure 2. The model shown in Figure 1 as 
applied to the magnetosphere. Depolarizing 
currents along high-latitude magnetic field 
lines close through the polar ionosphere. 
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In this view of the solar wind - magnetosphere 
interaction, geomagnetic activity is considered to 
be the magnetic effects of currents around and above 
the earth (.which may in turn induce currents inside 
the earth). These currents are partly associated 
with the depolarization of the magnetospheric bound- 
ary layer but to an even greater and more important 
extent they are induced by rapidly changing magnetic 
field configurations ( e.g . in magnetospheric fireball 
events) as the stressed magnetosphere gives way and 
relaxes to a lower energy state. As the basic inter- 
action is transfer of momentum from the solar wind 
to the magnetospheric and ionospheric plasma, we would 
quite generally expect that geomagnetic activity 
should increase with increasing momentum flux of the 
solar wind impinging on the earth's magnetic field. 

Any further solar wind parameters that could enhance 
the coupling to the magnetosphere might similarly be re- 
sponsible for additional enhancement of geomagnetic 
activity. 


2. RELEVANT PARAMETERS OF THE 
SOLAR WIND - MAGNETOSPHERE INTERACTION 

The initial entry of solar wind onto geomagnetic 
field lines seems to depend on the direction of the 
interplanetary magnetic field (Cole, 1974; Bahnsen 
and Hansen, 1976). In fact, the boundary layer is 
observed to be much thicker at times when the inter- 
planetary magnetic field has a large southward com- 
ponent (Sckopke et al. 3 1976). Over regions of the 
magnetopause where the magnetic fields outside and 
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inside the boundary are parallel, plasma flow is per- 
mitted across the boundary entailing an electric field 
tangential to the magnetopause. A plasma particle 
that has an initial guiding center velocity, carrying 
it across a boundary between magnetic fields of dif- 
ferent direction, will continue its motion when the 
two fields have a parallel component but will be re- 
flected back when the magnetic fields on either side 
of the boundary have an anti-parallel component. In 
loose terms we may say that the solar wind plasma can 
penetrate deeper into the geomagnetic field at places 
where the field direction is the same as the direction 
of the interplanetary 'magnetic field embedded in the 
solar wind because it takes longer for the plasma to 
realize that something is wrong. For anti-parallel 
fields, initial drifts towards the boundary actually 
result in a removal of plasma from the opposing mag- 
netic field lines. As pointed out by Cole (1974) 
this will tend to cause vacuum connection or merging 
of geomagnetic and interplanetary magnetic field lines. 
Vacuum connection is not associated with electric 
fields or particle energization. Turning a toy magnet 
in the earth's magnetic field causes a continuing 
change of the topology of the field around the toy 
magnet without electrical effects.. Solar wind plasma 
can now enter the magne-to sphere along the newly con- 
nected field lines as they are convected downstream 
and thus increase the thickness and extent of the 
boundary layer. There is experimental evidence that 
both the cross-field diffusion and the magnetic merg- 
ing discussed above are operating simultaneously 
(Reiff et al. 3 1977). In addition to the thickening 


7 



of the boundary layer, magnetic flux is also being 
transferred from the dayside magnetosphere to the 
tail region by the merging process. This in turn 
leads to increased frequency and intensity of re- 
connective events within the downstream magnetosphere 
resulting in enhanced geomagnetic activity. 

The number of interplanetary magnetic field lines 
that are brought up to the magnetosphere per unit 
time and unit area depends on the product of the mag- 
netic flux density B and the solar wind speed V, 

A geometrical factor depending on the angle between 
the interplanetary magnetic field and the geomagnet- 
ic field determines how much of the flux can merge. 

This factor is largest when the two fields are anti- 
parallel at the front of the magnetopause and decreases 
to zero for exactly parallel fields. Observations 
(Fairfield, 1967) show that the interplanetary mag- 
netic field is "draped" around the nose of the mag- 
netosphere as shown in Figure 3. The draping does 
not change the latitude angle of the field signifi- 
cantly. The result is then that the magnetic field 
just outside the magnetopause is tangential to the 
magnetosphere surface - directed from dawn to dusk 
in case of an ideal "away" polarity interplanetary 
magnetic field and from dusk to dawn in case of to- 
ward polarity. In addition the field may make a non- 
zero latitude angle with the ecliptic. We note that 
for negative latitude angles (southward fields) con- 
ditions are favorable for merging at the nose of the 
magnetosphere and also favorable for cross-field dif- 
fusion into the magnetosphere over the polar regions. 
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Figure 3- Equatorial plane view of the observed 
draping of the interplanetary magnetic field 
around the nose of the magnetosphere. The sit- 
uation is shown here for "away" polarity, 
the main effect is that the interplanetary mag- 
netic field at the dayside magnetosphere is 
largely parallel to the magnetopause. 


The geomagnetic dipole is roughly perpendicular to 
the solar wind velocity vector (implied in Figure 2). 
Qualitatively, it seems reasonable that the exclusion 
of solar wind plasma from the terrestrial environ- 
ment depends in some way on the angle il> between the 
solar wind flow direction and the dipole axis because 

the magnetic field seen by the solar wind is weakest 
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when ip = 9 0° . The angl-e ip varies both seasonally 
(90°+23?5) and diurnally (±11?5) due,to the 23?5 
angle between the earth's equator and the ecliptic 
and because the dipole axis is inclined 11? 5 to the 
rotation axis. It is thus likely that such varia- 
tions of the geometry of the magnetospheric "plasma 
physics experiment" could influence the .coupling 
efficiency. There are, in fact, observed variations 
in the amount of geomagnetic activity that closely 
follow variations of the dipole inclination ip . 

These are first of all' the classical semi-annual 
variation of the activity, first recognized by 
Sabine (1856), and* secondly the Universal Time var- 
iation fund by McIntosh (1959). Geomagnetic activity 
seems to be largest when the dipole axis is perpen- 
dicular to the solar wind flow direction. 

We have identified several solar wind - magneto- 
sphere parameters that are important for the genera- 
tion of and possibly modulation of geomagnetic activ- 
ity. In the following summary we discuss how funda- 
mental properties of the solar wind enter into the 
geomagnet ically active parameters. 

1) The density of the solar wind .momentum flux 
(the dynamic pressure) is given by pV 2 where p is 
the mass density and V is the solar wind speed. 
Often the proton number density, n, is used as a 
parameter largely proportional to p. The Helium 
content of the solar wind varies and can at times 
be quite high and contribute significantly to p. 

2) The influx of .merging interplanetary magnetic 
field lines is BV per unit length across the nose 
of the magnetosphere times a geometrical factor 
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q(a) where a is the angle between the directions 
of the draped interplanetary field and of the geo- 
magnetic field at the dayside magnetopause (see also 
Figure 18). The quantity B is the magnitude of 
the interplanetary magnetic field. The geometrical 
factor also includes the variation of the efficien- 
cy of cross-field diffusion of solar wind plasma 
with the angle a. It is worth pointing out that 
merging and cross-field diffusion takes place in 
different regions on the dayside magnetopause but 
that they both depend in approximately the same 
way on a. 

3) As geomagnetic activity and solar wind param- 
eters often are expressed as time averages over 
some interval, e.g. one hour or three hours, and 
since the relations between them are non-linear in 
many cases, the variability of the solar wind must 
be taken into account. This is particularly impor- 
tant for the direction of the interplanetary mag- 
netic field. It has been suggested that the varia- 
bility itself may contribute to geomagnetic activ- 
ity (Holzer and Reid, 1975). Because the time fluc- 
tuations are most pronounced for the direction of 
the interplanetary magnetic field, the following 
ratio will be taken as a measure of solar wind 
variability on the time scale involved (a few hours) 
f = o^/B , where 

a F = 0 X +ct Y +CT Z (3-) 

is the total variance computed as the sum of the 

variances for each component of the interplanetary 

magnetic field. 
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4) Finally we remark that the geometry of a dipo- 
lar magnetic field involves the ''colatitude" ^ in 
the form (1+3 cos 2 iji) which enters into the semi- 
empirical description of the variation of activity 
with the dipole inclination 

The ideal way of studying a phenomenon that depends 
on several parameters would be to let only one param- 
eter vary at. a time. This ideal can often be simu- 
lated if enough' observational data is available by 
grouping the data into several classes in such a way 
that within each class there is only a slight varia- 
tion of all the parameters. The number of classes may 
be decreased if some of the parameters are uncorre- 
lated (or only weakly correlated) in which case aver- 
age values of these parameters can be used. We shall 
employ the grouping technique extensively in the 
present study. 


3. GEOMAGNETIC ACTIVITY 

About one sunspot cycle worth of interplanetary solar 
wind data is available through the National Space 
Science Data Center (e.g., King, 1976). This includes 
both magnetic field data and plasma data-. To char- 
acterize geomagnetic activity, an index measuring the 
degree of disturbance is commonly used.. A great va- 
riety of geomagnetic indices have been proposed and 
employed over the years i some of them are specialized 
indices designed to characterize specific aspects of 
the total disturbance field while others are meant to 
be global or 'planetary' indices giving a measure' of 
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the worldwide level of magnetic disturbances (Siebert, 
1971; Mayaud, 1967 , 1968', 1972 ; Svalgaard, 1976), 

In the following section we discuss the meaning of 
the geomagnetic index am that will be used in the 
present invest iga'tion , Although the actual deriva- 
tion of such indices is a highly specialized subject, 
it is my experience that enough misconception and con- 
fusion exist about what geomagnetic activity is and 
how it is measured to warrant a somewhat tutorial 
approach (for complete details see Mayaud, 1967). 

Figure 4 shows a reproduction of a magnetogram of 
the variation of the horizontal component of the geo- 
magnetic field at a mid-latitude station. The record 
extends over twelve three-hour Universal Time intervals. 
A geomagnetic storm begins about 9 n on the first day 
and the field level the following night is strongly 
depressed relative to the level on the revious night. 

A regular daily variation is indicated by a dashed 
line as far as it can be determined from the undis- 
turbed portion of the record and from the observer's 
knowledge of what the typical daily variation gener- 
ally looks like for the station during the 'season (and 
lunar phase) in question. Winds and tides in the 
ionosphere are responsible for producing these rather 
regular - but nevertheless varying from day to day - 
excursions of the geomagnetic field. During each 
three-hour interval some irregular field variations 
appear in the record. They have an absolute maximum 
and an absolute minimum within the interval; part of 
this variation is due, however, to the regular daily 
variation and must be eliminated. Graphically, this 
corresponds to the measurement of the vertical dis- 
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Figure 4. Magnetogram of the horizontal component of the geomag- 
netic field at a mid-latitude- station. The record extends over 
twelve three-hour intervals. A magnetic storm begins around 9 h 
universal time. Dashed lines show the estimated daily variation. 
For several three-hour intervals the amplitude, a, of the activity 
is illustrated (after Mayaud, 1967). 




tance between two curves parallel ’to the daily vari- 
ation and enclosing the irregular variations. Several 
such amplitudes have been marked on Figure 4. It is 
important to emphasize that what is measured is not 
the deviation from some "quiet" reference level but an 
amplitude of fluctuations believed to be caused by the 
solar wind interacting with the magnetosphere during 
that three-hour interval. Finally, to construct a 
global or planetary index, amplitudes are averaged 
over selected observatories with as far as possible a 
uniform distribution in longitude in both the northern 
and the southern hemisphere. The details of the final 
computation (Siebert, 1971; Mayaud, 1968) need not con- 
cern us here; we only note that several indices exist 
because of differences in stations used and in time 
periods covered. Since 1959 the station distribution 
has been uniform enough to allow computation of a very 
close approximation to a truly global index: the am- 

index (Mayaud, 1968). the 'm' in am stands for mondial , 
the French word for worldwide. An earlier index, the 
well-known ap-index, (p = planetary) extends back to 
1932 but is based on an inferior station distribution. 
Using two antipodal observatories it has been possible 
to construct a homogeneous series of activity ampli- 
tudes going back to 1868: the aa-index (Mayaud, 1972). 

From the three-hour indices, daily indices can be 
computed. E.g., the Ap-index is the average ap-index 
over the UT-day. Expressing the Ap-index on a quasi- 
logarithmic scale from 0 to 9 gives a very compact re- 
presentation of geomagnetic activity on a daily basis: 
the C9-index. The- choice of a three-hour index as the 
basic activity measure is mainly dictated by the fre- 
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quency-spectrum of the irregular variations and may 
be justified on physical grounds because three hours 
is the time it takes the solar wind to pass the mag- 
netotail. The magnetosphere usually reacts to a 
change in boundary consitions in a time considerably 
smaller than three hours so that no time delay is nec- 
essary when comparing interplanetary data with the am- 
or ap - indices. 

4. INFLUENCE OF BULK PROPERTIES 
OF THE SOLAR WIND 

That geomagnetic activity indeed does depend on solar 
wind parameters, especially solar wind speed, was 
noted in the very earliest studies of the influence 
of the solar wind (Snyder et al. 3 1963), and has been 
extensively confirmed since. Figure 5 shows that on 
the average, the am-index is approximately proportion- 
al to the square of the solar wind speed V. It is 
somewhat remarkable that an index that was conceived 
almost forty .years ago (Bartels, 1938)' responds so 
clearly to an extra-terrestrial quantity. A vivid 
illustration of this relationship on a day-by-day basis 
is given by Sheeley et al. (1977). As we shall show, 
the relationship is really. a complex one; several phys- 
ical causes cooperate to give the overall effect. As 
the first step towards this goal we consider separately 
the possible influences of the influx of interplane- 
tary magnetic field lines, BV , and of the solar wind 
flux, nV 2 . Averages of the am-index computed for 
rather narrow bins or intervals of both BV and nV 2 are 
shown in Figure 6. We assume - and have actually veri- 
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am INDEX 



Figure 5. - Average dependence of the am- 
index on the solar wind speed V using data 
from the interval 1965-1973, More than 
9000 three-hour intervals of data are avail- 
able. The area of the filled circles are 
proportional to the number of three-hour 
intervals used in the average. For the 
larger datapoints the statistical uncer- 
tainty of the average is less than the dia- 
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meter of the circles. For an average over 
100 three-hour intervals the error is of the 
order of 10% of the average vaiues. A power- 
law dependence is suggested by the data with 
an exponent near 2.25. Similar results have 
been reported by Murayama and Hakamada (1975) 
and by Crooker and Feynman (1977). 


fied - that values of other solar wind parameters have 
approximately the same distribution in each bin. This 
is especially important in the case of the direction 
of interplanetary magnetic fields. Note that a doub- 
ling of BV also doubles the average am-index but that 
it is necessary to increase nV 2 by a factor of 7 or 8 
in order to obtain a doubling of am. The figure sug- 
gests the following approximate relationship 

am n. BV • (nV 2 ) x / 3 (2) 

implying that both the magnetic flux and the momentum 
flux are important -at the same time for the generation 
of geomagnetic activity. To eliminate the possibility 
that mutual correlations between V, B and n cooperate 
to produce this apparent relationship, we now examine 
only data with almost constant values of all param- 
eters except n and B. 

In the first- case we ask the question: If the num- 

ber density, n, is the only parameter allowed to vary, 
does the am-index vary as the one-third power of n? 

As shown in Figure 7, that seems indeed to be the case 
within ' observational accuracy. A similar analysis 
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Figure 6. Average values of the am-index for 
different ranges of the flux’ of interplanetary 
magnetic field lines, BV, and of the solar wind 
momentum flux nV 2 . The unit for BV in this 
figure is millivolt/meter. The number density 
n is given in protons/cm 3 and V q denotes the 
solar wind speed divided by 100: V Q = V/100, 

with V expressed in km/s. Each value on the 
plot represents an average of about 100 cases. 
Straight lines have been drawn corresponding to 
the relation given by eq.(2). 
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Figure 7. The dependence of am on the num- 
ber density, n, for almost constant values of 
all the other solar wind parameters. Filled 
circles represent cases where the interplan- 
etary magnetic field was largely southwards 
with respect to the geomagnetic field (i.e. 
cosa<-0.25); open circles represent cases of 
predominantly northwards fields (cosct>+0 . 2'5) , 
while circles with a dot represent cases of 
interplanetary fields largely perpendicular 
to the geomagnetic dipole. A typical error 
bar is shown in the center of the figure. 
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n =5.0+ 1.0 cm 
V=395± 55 km/S 
cr/B = 0,5± 0.1 
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Figure 8. Similar to Figure 7 but in this 
case the field magnitude B is used as the 
independent variable. In both Figures 7 and 
8 the selected ranges of solar wind param- 
eters are shown in the upper left-hand corner 



for the magnetic field strength, B, is shown in Fig- 
ure 8, confirming the first-power dependence. In both 
figures a further subdivision according to the angle, 
a, between the direction of the geomagnetic field on 
the dayside magnetosphere and of the interplanetary 
magnetic field was performed. The data seems to be 
consistent with a first-power dependence on B and a 
one-third power dependence on the number density n. 

The relationship eq.(2) indicates that the influences 
of those two parameters are not independent -- in . 
that case the effects would be additive rather than 

i 

multiplicative. In addition it is apparent that the 
direction of the interplanetary magnetic field is of 
paramount importance and that its influence can be 
described as a modulating angle-dependent factor that 
should be incorporated into eq.(2) - because of the 
roughly parallel sequences of datapoints in both Fig- 
ures 7 and 8. 

That geomagnetic activity correlates very well with 
BV times a function q(a) of the angle a between the 
interplanetary and the geomagnetic field is well-known 
and generally accepted (e.g., .Burton et al. 3 1975 ; 
Crooker, 1975). In fact, one of the outstanding cur- 
rent problems in- the study of the magnetosphere is to 
understand- the physics, of that correlation. By form- 
ing the ratio am/BV, we get a quantity that is re- 
duced for the influence of BV and thus might show a 
clearer correlation with the momentum flux nV 2 . The 
result is shown in Figure 9. A physical interpre- 
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tation of this relationship might be that the magneto 
sphere is powered by solar wind kinetic energy, but 
that the power input is strongly controlled by cou- 
pling via the interplanetary magnetic field. 




am’ 


Figure 9. By dividing the am-index by BV we 
isolate the influence due to nV 2 . In this 
figure, am/BV o is plotted against nV 2 (where 
V q = V/100) confirming the one-third power 
dependence suggested in Figure 7. The areas 
of the filled circles are proportional to 
the number of data values in each average. 

A total of more than 9000 three-hour intervals 
were available for the analysis. 
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5. INFLUENCE .OF GEOMETRICAL FACTORS 

Further insight into the coupling problem can be gain- 
ed by performing the same type of analysis as in Fig- 
ure 9 for various subsets of the data characterized by 
different ranges of the relevant parameters such as 
the angle a between the two fields, the relative vari- 
ability f and the angle tf> between the solar wind 
flow and the earth's dipole axis. Figures 10, 11 and 
12 display the results. The effect of varying the 
angle a is by far the most important (Figure 10), al- 
though significant increase in activity is also asso- 
ciated with high variability of the interplanetary 
magnetic field (Figure 11). Finally there is a small 
but persistent tendency for activity to be higher when 
the geomagnetic dipole is perpendicular to the solar 
wind flow direction (Figure 12), The consistent 
trends in these Figures strongly suggests that we may 
largely remove the influence of BV and nV 2 simply by 
computing a reduced am-index: 


„ t _ /<BV>\ 

am' = am (— ■ ) 

H >) 1/3 

(3) 

\ BV / 

\ nV 2 / 



We have used <BV> = 2100 and <nV 2 > = 1.05 x 10 6 (B in 
nT, V in km/s , and n in protons/cm ). To investigate 
the effect of a and of the relative variability f in 
detail, the data is divided into groups according to 
the value of f. The group characterized by f = 0.5 
will consist of all data values for which f lies be- 
tween 0.45 and 0.55, and so forth. Then for each 
group we make a further subdivision according to the 
value of cos a . There is enough data to allow a bin- 
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width of 0.1 in cosa without compromising the accu- 
racy of the mean value of the am-index computed for 
every subdivision. Figures 13 and 14 show the result. 



Figure 10. Similar in format to Figure 9, ex- 
cept with a further subdivision of the data into 
five classes of different orientation of the in- 
terplanetary magnetic field; from mostly south- 
wards (filled circles) turning to mostly north- 
wards (open circles). The class intervals se- 
lected for cosa are shown in the upper left cor- 
ner . 
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Figure 11. Similar to Figure 10, except 
that the data was divided according to the 
relative variability f = a/B of the inter- 
planetary magnetic field. The classes 
range from quiet (f<0.37; open circles) to 
disturbed {f>0.60,; filled circles) con- 
ditions. Note that the time scale f.or 
determination of f is three hours. If 
the field magnitude iwas constant {very 
nearly -true) within a three-hour interval 
the maximum value of f >would be 1. 
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25 50 100 200 300 400 

nV 0 2 

Figure 12. Similar to Figure 11, except 
that the- data was divided in two classes 
according to whether the tilt angle (y = 90°-i|i) 
of the geomagnetic dipole was numerically 
less than 14° (filled circles) or greater 
than 14° (open circles). There is clearly 
a persistent tendency for activity to be 
higher when y is small, i.e. when the di- 
pole is nearly perpendicular to the solar 
wind flow direction. 
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Figure 13. The dependence of geomagnetic 
activity Cam-index reduced for BV and nV 2 
influence, see text) on cosa, where a is the 
angle between the geomagnetic field lines at 
the nose of the magnetosphere and the aver- 
age interplanetary magnetic field during 
the three-hour interval corresponding to am. 
Draping (Figure 3) of the interplanetary 
field is assumed and the reader is referred 
to Figure 18 for a precise definition- of a. 
Several curves are shown for various values 
of the relative variability f; from f = 0.0 
to f = 0‘. 9 as labeled on the- figure. Note 
the very flat curves for high variability 
and the steeply falling curves for low vari- 
ability of the interplanetary field. 
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Figure 14. Dependence of reduced am-index 
on the relative variability f of the inter- - 
planetary magnetic field for several val- 
ues of cosine of the angle a between the 
interplanetary field and the geomagnetic 
field. Note that for most values of a ac- 
tivity increases with increasing variability, 
but that for nearly southwards fields the 
inverse is the case. The curve labeled 0 
is for cosa = 0 which is the most probable 
value of cosa and is thus indicative of 
the typical situation. Figure 16 shows 
some of the individual data points used to 
construct the smooth curves in Figures 13 
and 14. 
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arm= am ( f, cos a*) 



Figure IS. Actually Figures 13 and 14 are - 
synthesized from this figure. Bins of cosa, 
(0.1 wide) and bins of f (Oil wide) were form-r 
ed and the average value of reduced am was 
computed for each bin. The figure shows a 
contouplot of the results. Contour levels 
less than average are shown as dashed lines. 
It is now a simple matter to construct the 
curves in Figures 13 and 14. The following 
fourth-order polynomial fit to the contours 
gives synthetic am-values that have a very 
high correlation (0.992) with the observed 
averages for each bin: 
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am = 6.6 exp ( -1 . 275cosct + l . 815/+0 . 565cos 2 a+ 

+ 2 . 0 41/cos a- 0 . 642/ 2 -0 . 001cos 3 a-l . 877/cos 2 a- 
-2.715/ 2 cosa + 0. 63 6/ 3 + 0.118cos 4 a+0. 0 81/cos 3 ct + 
+1 . 190/ 2 cos 2 a+l . 93 5 /cos a- 0 . 754/ 4 ) 

The above formula expresses a purely formal 
relationship, of course. 



cos a 


Figure 16. Reduced am as a function of cosa 
for f = 0.3 (open circles), f = 0.4 (circles 
open on the left), f = 0.5 (circles filled on 
the left), and f = 0.6 (filled circles). The 
scale for am is logarithmic. An expression of 

the form am = kexp(-gcosa) fits the data. 
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Each curve in Figure 13 is drawn through the average 
points for each f-group and shows how the activity 
varies with a for that particular variability of the 
interplanetary magnetic field. In Figure 14 we show 
how the activity varies with the variability f for 
different values of cosa ranging from -1 to +1, i.e. 
from southwards to northwards. 

We note first that when the variability is high 
(e.g. a f = 0.9) geomagnetic activity is almost inde- 
pendent of a. The reason is, of course, that the av- 
erage direction of the interplanetary magnetic field 
during the three-hour interval in question is unde- 
fined or nearly so. With decreasing variance the 
importance of southwards fields (cosa<0) becomes 
more and more prominent. For f = 0 the curve rises 
sharply as cosa approaches -1. In discussing the 
physics of this relation we should bear in mind that 
the interplanetary magnetic field at the magneto- 
pause fluctuates more than the field in interplane- 
tary space due to passage through the bow-shock of 
the magnetosphere. This suggests that the curve 
showing am-index dependence on cosa extrapolated to 
zero variability of the field at the magnetopause 
would rise even sharper than the curve labeled f=0 
in Figure 13. The implication seems to be that the 
coupling to the solar wind due to magnetic field con- 
nection is very weak unless the geometry is very fa- 
vorable, i.e. the external field is almost anti-par- 
allel to the dayside geomagnetic field. Due to ever- 
present fluctuations of the interplanetary magnetic 
field - considerably enhanced after passage through 
the bow-shock - favorable conditions for connection 
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occur often enough within a three-hour interval at so 
many places on the magnetopause as to give the impres- 
sion that reconnection and hence geomagnetic activity 
occur for all orientations of the ' interplanetary mag- 
netic field and varying in efficiency smoothly from a 
maximum for anti-parallel fields to a non-vanishing 
minimum for parallel fields. 

A simple computer simulation can qualitatively re- 
produce the curves in Figures 13 and 14 starting with 
the f = 0 curve and by just varying the three compo- 
nents of the interplanetary field such as to result in 
a certain average vector making the angle a with the 
geomagnetic field. By adjusting the time scale of the 
computer generated variations, any relative variabili- 
ty f = a/B can be produced. Finally by averaging am- 
values taken from the f = 0 curve corresponding to 
each of the "f ine-scale" values of a, synthetic curves 
can be constructed showing the average am-index for 
given cosa and given f. The results are qualitatively 
very similar to what is shown in Figures 13 through 15 
By assuming that the variability of the north-south 
component is considerably greater than that of the 
other components of the magnetosheath field one can 
even obtain reasonable quantitative agreement. 

As our knowledge of fluctuations of the field in 
the magnetosheath - that turbulent layer between the 
bow-shock and the magnetopause - is rather limited 
(e.g. 3 Fairfield, 1976) it is not possible at the pres 
ent to extend the simulations mentioned above to more 
realistic conditions, but it already seems likely that 
most - if not all - of the difference between the de- 
pendence of am on cosa for different degrees of varia- 
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bility can be understood simply in terms of aver- 
aging the non-linear response of am to cosa over 
varying amounts of variance. It does not seem neces- 
sary to invoke any ideas that the variance per se is 
instrumental in producing geomagnetic activity. As a 
matter of fact, as is evident from Figure 14, for 
southwards fields, activity actually decreases with 
increasing variability of the field. 

The relative variability, f, is typically near 0.5. 

In Figure 16 we have plotted the value of reduced am- 
index as a function of cosa for several values of f 
near 0.5. The scale for am is logarithmic and the 
resulting nearly linear trends of the datapoints sug- 
gest a simple exponential relation 

am = k(f )e~ gCf )cosa (4) 

where k and g vary slightly with f. Within the typical 
range of f we may consider g to approximately constant, 
being near 3/4. Although such an empirical relation 
is very convenient to work with, we should emphasize 
that it is just that and that it not immediately im- 
plies physical reality or insight. The value of k fol- 
lows from Figure 14 by setting cosa equal to zero. 

6. COMPUTATION OF THE AM-INDEX 
FROM SOLAR WIND DATA 

Statistically, the variability, the field magnitude, 
and the plasma density all depend slightly on the so- 
lar wind speed. Figure 17 shows the average relative 
variability as a 'function of V. A similar increase 
with V is found on the average for the field strength, 
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V (km/s) 


Figure 17. Average value of the relative vari- 
ability f = cr/B of the interplanetary magnetic 
field as a function of solar wind speed V, com- 
puted for 25 km/s bins. The histogram shows the 
number of three-hour intervals in each bin. 

A linear relation f = 0.28+ 4.2*10~ 4 V fits the 
data. Over the observed range of V, a power-law: 
f^v 0 ' 4 is an equally good fit. Dividing the data 
into two parts: dV/dt>0 and dV/dt<0 gives essen- 
tially the same result for each partition except 
that f is 10% higher for the subset with increas- 
ing velocity than for the subset with decreasing 
speed. A similar result holds for the field mag- 
nitude B although somewhat more noisy. 
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p <0 


p >0 




AWAY 


TOWARD 


. °A = 


a T = 180- A-/? 

- = 180-{A+/3) 


Figure 18. Field line geometry at the nose of 
the magnetosphere as seen from the sun for dif- 
ferent orientations of the interplanetary mag- 
netic field (dashed arrow) and the geomagnetic 
field (solid arrpw marked SN). The ecliptic is 
indicated by EE. The interplanetary field is 
assumed to be ’'draped" around .the nose of the 
magnetosphere as shown in Figure 3. This means 
that away-polar ity is directed from dawn to dusk 
(left panel) and toward polarity is directed from 
dusk to dawn (right panel). The angle A is the 
angle between the geomagnetic field lines and the 
ecliptic measured from the ecliptic towards the 
field lines and varying from 55 to 125 depending 
on season and on universal time. 8 is the lati- 
tude angle of the interplanetary field. The 
angle a between the two fields is now defined as 
shown separately for the two polarities. 
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while the density is anti-correlated with V. For the 
range of V from 250 km/s to 850 km/s all these rela- 
tionships can be approximated with power-laws: 

f^V 0 * 4 , k^V 0,25 , B'W °* 4 , no-V -0 * 5 (5) 

Combining (2) through (5) under typical conditions 
(cosa=0) we get that the am-index should 'depend on 
solar wind speed approximately as 

am ^ BVCnV 2 ) ^ * k V 2,15 (6) 

We may compare this result with Figure 5 and note rea- 
sonable agreement. 

We have thus found that the rather simple state- 
ment that geomagnetic activity increases roughly with 
the square of. the solar wind speed does not reflect an 
equally simple physical situation but is a consequence 
of several physical and statistical conditions cooper- 
ating in a rather complex manner. The major contribu- 
tions arise from the BV-factor and from the momentum . 
flux factor (nV 2 ) to the one-third power. The physics 
of the specific way the momentum flux or the dynamic . 
pressure of the solar wind enters into the functional 
expression for geomagnetic activity such as eq.(2) is 
not well understood. The same is the case with the 
details of magnetic reconnection or plasma cross- 
field diffusion. On the other hand, we have obtained 
a rather detained understanding - or at least a 
quantative description - of how solar wind parameters 
affect the geomagnetic field on the larger scale. 

The degree to which we can reproduce the observed am- 
index from solar wind data may be judged from Figure 
19 which shows for a number of Bartels rotations both 
observed and reconstructed am-indices. The comparison 
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OBSERVED AND COMPUTED AM-1ND1CES 



3 6 9 12 15 18 21 24 27 

DAY IN SOLAR ROTATION 


Figure 19 


38 






Figure 19 . Observed and reconstructed am-indices 
for Bartels rotations 1836-1841 (Oct. 3, 1967 to 
March 12, 1968). For every three-hour interval 
where solar wind data was available, am was com- 
puted using 


am 


6.6 q(f ,a) 


JBV 0 

(21 


nVp 2 ] 1 / 3 , 1.157 

- (1+3cos 2 iJi)2/3 


where V 0 = V/100 (V is solar wind speed in km/s), B is 
field magnitude in nT, and n is the number den- 
sity in protons per cm 3 . The function q(f,a) of 
relative variability f and angle a is the fit- to 
reduced am given in the caption to Figure 15. 

The scale of am is logarithmic because we want to 
verify the synthesized am-indices against obser- 
vations over the full range of the index rather 
than just for the higher values of the index. 

The two overlapping curves show the two .indices 
for times when solar wind data was available. 

Where only a single curve is visible over a cer- 
tain time interval it just means that the computed 
and the observed indices track each other perfectly. 
The main area of disagreement is for very small 
am-values (am<l; all such cases are plotted as 
am=l) . 


Figure 20. a) Computed values of cos(A-6 + ) for 
each three-hour interval of every month. -The re- 
sult is presented as a contourplot (negative con- 
tours are dashed) . The angle A depends on day of 
year and on Universal Time. Let X be the ecliptic 
longitude of the mean sun and let d denote day of 
year (Jan. 1 h 1) . We then have the following ap- 
proximate relations for computation of A: 
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(Figure 20, continued) 


A(d) = 279?69 + 360°/365.24 (d-h) 
?z (UT) = 360°/24 (UT-10' h 40 m ) + X 


cosi|/ = sinAcosisine - sin-£(sinAcosecos7i - cosAsin h) 
cosx = cosAcos-£sine ** sinf (cosAcosecosfc + sinAsinft) 

and finally 


cosA = cosx/sini^ 

e = 23?5 is the obliquity of the ecliptic and 
i - 11? 5 is the geographical co-latitude of the 
earth's magnetic pole. 

b) Observed yearly and Universal- Time variation 

of ^ 

S(t) = 2 • (am, - am m j)/(am. + am m ,) 

Away Toward Away Toward 

for the interval 1963-1974 when interplanetary mag- 
netic field data was available. For each month 
and every three-hour UT-interval we find the aver- 
age values of the am-index separately for Away po- 
larity and for Toward polarity. More precisely, 
the sign of the azimuthal component of the inter- 
planetary magnetic field has been used as the se- 
lection criterion (By-O = Away; By<0 = Towards) be- 
cause it is this component that controls the way 
the interplanetary magnetic field is draped around 
the nose of the magnetosphere (see Figure 3). 

A slight smoothing of the data has been performed 
by the contour plotting routine. 
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is quite satisfactory over the full range of activity 
from quiet to very disturbed. Similar comparisons have 
been carried out for the entire sunspot cycle with -the 
same- -sa-ti-s-factory result . The main' area of disagree- 
ment seems to be f,or extremely quiet intervals (e,g. f 
am = .0) where the reconstructed values turn out a 
little too high - typically about 2 rather than 0. 
Possibly the curves for coset near +1 in Figure 14 
should be. lowered somewhat or maybe the am-index is a 
poorer measure of geomagnetie activity for very small 
values of the index where it is extremely difficult 
to measure . 

7. SEMI-ANNUAL AND UNIVERSAL TIME VARIATIONS 

The usefulness of the concept of the earth being a 
probe for interplanetary and solar conditions depends 
on our understanding of the "instrumental response 
function!’. The foregoing analysis suggests -that we 
have reached a point where physically meaningful 
quantitative results can be extracted from geomag- 
netic data. Let us consider an idealized interplan- 
etary magnetic field configuration. The field is 
wound by solar rotation into a spiral in the equato- 
rial plane of the sun. Because the ecliptic is in- 
clined 7%°to the solar equator the ecliptic latitude 
angle of the average interplanetary magnetic field 
measured at the earth vary with time of year: 

3 = ~P 3 0 cos t (7) 

where 3 0 = 7?25 and t is reckoned from 0 to 2ir from 
December 7 when the heliographic latitude of the 
earth is zero. Note that for the two opposite po- 
larities of the field g has opposite signs (p=+l = 
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polarity away from the sun; p=-l = polarity towards 
the sun) . If A is the angle between the ecliptic and 
the geomagnetic field lines at the nose of the magne- 
tosphere (see Figure 18), then 

<x + = A-B + , ct_ = 180°-(A+Bj (8) 


thus 

cosa = -eos(A+8_) = -cos(A~8 + ) = -cosa + (9) 


Using eq.(4) we may now compute the expected average 
value, a, of geomagnetic activity assuming equal prob- 
ability of both polarities: 

-gcosa + +gcosa + 
a(t) = %(a + +a_) = ^kCe + e ) 

a(t) = k*cosh (gcos(A-8 + ) ) (10) 

Similarly, we get for the difference in activity be- 
tween the two polarities 


Aa(t) = a + -a_ = -2k*sinh(gcos(A-B + ) ) (11) 


The ratio 6(t) = Aa(t)/a(t) is independent of k and 
depends only on the geometry of the situation 


5 ( t ) = -2 *tanh(gcos(A-3 + ) ) (12) 

Because gcos(A-8 + ) is always numerically less than 0.5 
we may write 

5 ( t ) = -2g*cos(A-8 + ) (13) 

In Figure 20a we show the variation of cos(A-6 + ) with 
time t, computed for. each three-hour UT-interval of 
every month of the year; and in Figure 20b we show for 
comparison the observed variation of fi(t). The agree- 
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ment is quite good with respect to phase; the amplitudes 
of the variations shown in the figure are very nearly, 
the' same, corresponding to a value of 2g of unity, or 
g = 0/5. 'This is somewhat smaller than the 0.75 found 
earlier using individual three-hourly observations 
rather than assuming an' ideal field configuration. In 
fact we get here a measure of the error committed by 
working with the idealized field instead of the actual 
field orientation. 

The Universal Time variation of 6(t) has a maximum 
at 10^30 m UT and a minimum at 22 il 30 m ' i while the yearly 
variation has a minimum in early April and a maximum in 
early October.' The amplitudes of these changes are 
very substantial amounting to 60% attesting to the 
strong control of the activity by the field orientation. 
Because cosh(x) = cosh(-x), eq.(10) predicts a semi- 
annual variation of a(t) as was first realized by 
Russell and McPherron (1973), but this variation is 
of insignificant amplitude and is almost completely 
masked by another semiannual variation related to the 
dipole tilt angle ( e.g. 3 Berthelier, 1976), Figure 21 
is contourplots of observed average am-index in the 
same format as used in Figure 20. Furthermore, high 
activity is indicated by dotted regions and is pre- 
dominantly observed near the equinoxes. The amplitude 
of the observed variation is three times larger than 
that predicted by eq.(10) and the phase is very dif- 
ferent. Eq . (10 ) predicts maxima at times when 5(t) 
attains extremum values ( c.f . Figure 20) quite contrary 
to the observed phase of the average am-variatioh ( c.f . 
Mayaud, 1974a). We are then forced to concede that a 
further modulation of geomagnetic activity exists pro- 
ducing the variations sho.wn in Figure 2lb and the dif- 
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ference between the two lines in Figure 12. As first 
suggested by Bartels (1925) it seems that geomagnetic 
activity is enhanced when the earth’s dipole axis is * 
perpendicular to the solar wind direction. That we are 
dealing with a modulation rather than an excitation 
mechanism is suggested by the basically parallel lines 
in Figure 12 showing that the influence of the dipole 
tilt is proportional to the level of activity itself. 


Figure 21. a) Contourplot computed from eq.(16) 
with am 0 set equal to the average value of am. 
over the 1959-1974 period (see below), 
b) Contourplot of the observed variation of the 
am-index with time of year and with Universal 
Time. The observed variation is based on 16 years 
of data (1959-1974). Regions of maxima are dot- 
ted. The first full contourline corresponds to 
am = 21, and lines are drawn one am-unit apart. 
Contour lines in regions of lower than average 
values are shown as dashed lines. 

The correlation between the two contour plots is 
0.90, i.e. a) is a very good fit to b) . 

In this figure no separation according to polarity 
has been performed so we should have a plot of the 
quantity a(t), i.e. activity averaged over both 
polarities. Note that the UT-variations for the 
two solstices are in anti-phase and that the UT- 
variation at the equinoxes i’s small with two maxima 
and two minima during the UT-day. Other magnetic 
indices ( e.g.j Kp or AE) show similar yearly varia- 
tions but very distorted UT-variations due to un- 
even station distributions. 
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Figure ?1 
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The physical explanation for the dipole tilt mod- 
ulation of activity is still obscure. It has been 
suggested (Boiler and Stolov, 1970; Mayaud , 1977) 
that the Kelvin-Helmholtz instability of the boundary 
between two magnetohydrodynamic fluids in relative 
motion may be operating at the flanks of the magneto- 
sphere leading to release of some of the energy 
stored in the magnetotail. An approximate instability 
criterion may be expressed as 


M 2 > 1 + (B./B ) 2 cos 2 if (14) 

A i o 

where B^ and B^ are the magnetic field strengths just 
outside and just i«side the magnetopause. If the 
Alfvenic Mach-number = (solar wind speed) / (Alfvdn 
speed) exceeds some value depending on the angle 
between the solar wind direction and the dipole axis 
then the magnetopause may be unstable against the 
growth of the Kelvin-Helmholtz. instability, thereby 
maybe resulting in a greater stress on the magneto- 
sphere . 

It is also possible that the coupling between the 
solar wind and the magnetosphere is dependent on the 
basic geometry of the magnetic field around a dipole: 

B = B (1 + 3cos 2 ^)^ (15) 

• eq 

where is the field strength in the equatorial plane. 

In any case, the parameter (l+3cos 2 ifO has been found 
to enter into a simple empirical expression for the 
modulation, namely 


am 


1.157 am D (l + 3cos 2 i|0 


-2/3 


(16) 


which produces an excellent fit to the observed am- 
variation with Figure 21a shows the variation of 
the above relation with time of year and time of UT- 
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day. The constant 1.157 is just a .normalization fac- 
tor (-equal to the yearly mean value of (l+3cos 2 $) 2/|f3 ) 
and am 0 is the unmodulated value of am. Of course, the 
relation given by eq.(16) is not unique, but there can 
be little doubt that some' function of accounts for 
the observed variations in considerable quantitative 
detail. This result puts an upper limit on systematic 
variation of solar wind speed with heliographic lati- 
tude. Judging from Figure 21, it seems that no more 
than two, am-units (nT) are not already accounted for. 
This corresponds to about 10%. Referring to Figure 5, 
we can translate that into a maximum variation with 
heliographic latitude of the wind speed of 16 km/s over 
7 degrees of latitude, or 2,3 km/s per degree. Any 
variation larger than that would distort the variation 
shown in Figure 21 in an unacceptable way. 

8. VARIATIONS WITH HELIOGRAPHIC LATITUDE 

In a careful analysis of a 100-year series of the geo- 
magnetic activity index aa, Mayaud (1974b) showed that 
the semi-annual variation attains maxima 4. 0+2. 5 days 
after the equinoxes. This difference can be understood 
in terms of the aberration of the solar wind caused by 
the earth's orbital movement. A significant variation 
of solar wind speed with heliographic latitude would 
shift the time of maxima towards March 7 and September 
7 when the earth is at the highest latitudes. No such 
shifts are observed when activity series of sufficient 
length are analyzed. If high-speed solar wind streams 
predominantly originate in .coronal holes that are lower- 
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latitude extensions of the polar holes ("polar cap lobes") 
one would naturally expect at least a statistical tendency 
for solar wind speed to be higher at higher heliographic 
latitudes. The situation is analogous to the - classical 
axial-equinoctial controversy among students of terrestrial 
magnetism. The axial explanation of the semiannual vari- 
ation takes its starting point in the fact that sunspots 
occur most frequently in heliographic latitudes higher 
than 10° while the equatorial belt of the sun is rela- 
tively free of spots. If, then, the solar wind streams 
originate in the same belts in which the spots occur 
they should more likely sweep across the earth in March 
and September causing an enhancement of geomagnetic 
activity at those times of the year. This explanation 
is very simple and intuitively appealing in contrast 
to the equinoctial hypothesis that ascribes the en- 
hancement to an obscure and largely unspecified mech- 
anism whose efficiency is somehow - in an unknown way- 
modulated by the angle between the geomagnetic dipole 
and the solar wind flow direction. Today - as forty- 
five years ago (Bartels, 1932) - the axial hypothesis 
must be discarded because of its discordance with obser- 
vations. Direct in~situ observations of solar wind 
speed (Bame et at. ., 1977) and of interplanetary magnetic 
field strength (King, 1976) over one full sunspot cycle 
show indeed no detectable variation of either parameter 
with heliographic latitude within ±7° of the solar 
equator . 

Such a variation is, however, found in the polarity 
of the interplanetary magnetic field. Rosenberg and 
Coleman (1969) found that the probability of occurrence 
of a given polarity in a 27-day solar rotation varies 
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in a roughly sinusoidal way during the year in phase 
with the heliographic latitude of the observer (at or 
near the earth). Wilcox and Scherrer (1972) extended 
the analysis back in time to cover more than four sun- 
spot cycles and showed that the phase of the above 
’Rosenberg-Coleman * effect reverses 2-3 years after 
each sunspot maximum coinciding with the inferred 
times of reversals of the solar polar fields. The 
Rosenberg-Coleman effect owes its existence to the 
fact that the interplanetary magnetic field within 
several astronomical units of the sun appears to have 
one polarity in most of the hemisphere north of the 
solar equatorial plane and the opposite polarity in 
most of the hemisphere south of the equatorial plane 
(Svalgaard and Wilcox, 1976; Smith et al. 3 1977). The 
two hemispheres are separated by a warped or curved 
current sheet that typically crosses the solar equa- 
torial plane in 2 to 4 places, thus dividing the equa- 
torial region into 2 to 4 sectors of alternating polarity. 
Near sunspot minimum the latitudinal extent of the warps 
of the curved current sheet is of the order of 10°, so 
that the sector boundary (f.e., the current sheet sepa- 
rating the two hemispheres of opposed field polarity) 
is almost parallel to the solar equatorial plane. It 
is then clear that at such times going out of the solar 
equatorial plane will increase the probabilty of obser- 
ving a .certain polarity over that of observing the op- 
posite polarity. If the observer reaches a latitude 
greater than the extents of the sector-warps in the 
current sheet, he will see a unipolar interplanetary 
magnetic field with no sector structure at all (Smith, 
et al . , 1977). 
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9. THE 2 2 -YEAR CYCLE IN GEOMAGNETIC ACTIVITY 


Let us model the Rosenberg-Coleman effect with an ex- 
pression of the form: 

w - h + p»r*sin t (17) 

P 

where w^ is the probability of observing the polarity p. 
Time t varies from 0 to 2 tt during the course of a year 
and is reckoned from December 7 when the heliographic 
latitude of the earth is zero. The magnitude of the 
coefficient r depends essentially on the latitudinal 
extent of the warps in the sector boundary current 
sheet, while the sign of r changes when the polarity 
of the solar poles reverse shortly after each .sunspot 
maximum. Figure 22 shows the value of r as a function 
of the warping of boundary. 

Geomagnetic activity averaged over both polarities 
may now be written 

a(t) = w + a + + w_a = %(a + +a_) + r(a + -a_)sin t 

a(t) = a(t) (l+r<5 (t)sin t) (18) 

In order to simplify the considerations we will work 
with daily averages so that the Universal Time varia- 
tion of 6 ( t ) is averaged out. In that case we find 
that 6 ( t ) varies very nearly sinusoidal during the 
year but with a phase that is different from the phase 
of the heliographic latitude variation. Let the phase 
difference be n (n~33 > ) then from eq.(13) and Figure 20: 

6(t) = -2g* 0 .43sin(t+n) (19) 

Averaged over a full year sin( t ) sin( t+n ) is equal to 
^scosri so that the yearly average of a(t) can be written 
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LATITUDINAL EXTENT OF SECTOR BOUNDARY WARPS 

Figure 22. Amplitude, r, of the Rosenberg- 
Coleman effect at the earth as a function of 
latitudinal extent of the curved current sheet 
that separates opposite field polarities in in- 
terplanetary space (see text for definition of 
r). The number of sectors per rotation (2, 4, 
or 6) has a slight influence on r as shown. 
Computation of r was done here by a simple com- 
puter simulation of the geometrical situation. 
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GEOMAGNETIC ACTIVITY INDEX aa 



YEAR OF SUNSPOT CYCLE 


Figure 23. Variation of yearly means of geomag- 
netic activity as a function of phase within the 
sunspot cycle. Activity is measured by the aa- 
index which for this purpose can be considered 
identical to the am-index. Even-numbered (filled 
circles) and odd-numbered (open circles) cycles 
are considered separately. At the bottom is shown 
the average sunspot numbers for the solar cycles 
since 1868 when the aa-series begins. Note that 
during the first half of the cycle, odd cycles are 
more geomagnetically active than even cycles while 
the reverse is true during the last half of a sun- 
spot cycle. 
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( 20 ) 


<a> = <a>(l-0.43gr*cosn) = <a->(l-0 . 28r ) 

year - ■ 

Thus for the eleven years. - .f.r-om short-l-y after one. sun- 
spot maximum to shortly after the next - when r is neg- 
ative geomagnetic activity will be increased, while for 
the next' eleven years when r is positive activity will 
.be diminished. This 22-year cycle in geomagnetic act- 
ivity is actually observed (Chernosky, 1966; Russell; 
1974) as shown in Figure 23, depicting the average sun- 
spot cycle variation of geomagnetic activity shown sep- 
arately for odd- and even -number e.d sunspot cycles. 

If we set r = 0.25, corresponding to a maximum sector 
■boundary warp of 10°, the difference in activity between 
odd and even cycles should be about 14% according to 
eq.(20) to be compared with the 20% actually observed. 
Considering the many simplifications we have intro- 
duced, the agreement is quite satisfactory. The 22-year 
cycle is a good example of how geomagnetic measure- 
ments can- be a valuable supplement to direct obser- 
vations of the sun — polar field reversals have only 
been directly observed twice. 

10. VARIATIONS OF GEOMAGNETIC 
DISTURBANCES WITH SOLAR ACTIVITY ’ 

When it was thought that geomagnetic activity was caused 
by corpuscular streams emitted from active regions on 
the sun (either sunspots or facular areas), the sunspot 
cycle variation of geomagnetic activity was taken for 
granted, as something that hardly required any elabo- 
ration. Today, the primary effect of solar activity is 
generally believed to be caused by compression of the 
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ambient solar wind by blastwaves emitted from solar 
flares. The ensuing high values of solar wind speed, 
density and especially magnetic field strength ensure 
major events of violent activity: geomagnetic storms 

which have an obvious sunspot cycle dependence. But, 
in addition, variations of fundamental solar wind para- 
meters can be very important. The 1973-75 period re- 
sulted in exceptionally high geomagnetic activity to 
the extent that it is hard to speak about a- regular 
sunspot cycle variation of geomagnetic activity for the 
period 1965-1976. In a similar way solar wind speed 
has had a pronounced tendency to have higher values 
in the years preceding sunspot minimum ever since space- 
craft data became available in the early 1960’s 
(Gosling et al. 3 1976; Bame et at ., 1977). On the other 
hand, no sunspot cycle related variation of magnetic 
field magnitude has been found (King, 1976). In order 
to put the recent data in proper perspective it Is 
necessary to consult the historical record. Geomag- 
netic activity has been monitored for almost two 
centuries and reliable indices of constant calibration 
exist for little longer than the past one hundred years. 
As shown in Figure 24, very significant long-term trends 
exist in the activity record. The sunspot cycle vari- 
ation is generally discernible in addition to a vari- 
ation on a time scale of 80-100 years. It is inter- 
esting that the last sunspot cycle is not at all typical 
(although maybe #16 from 1924 to 1933 is somewhat sim- 
ilar to cycle #20). We should, therefore, exercise 
caution in drawing general conclusions about the sun- 
spot cycle behavior of the solar wind parameters from 
observations covering only the last cycle. 
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Figure 24. Yearly means of geomagnetic activity 
from 1781 to the present. From 1868 the index is 
homogeneous and has constant calibration as dis- 
cussed by Mayaud (1972). The earlier data has been 
derived from the daily ranges of magnetic declina- 
tion as given by Wolf (1884). This earlier series 
covers the interval 1781 through 1880. For the in- 
terval 1868-1880 the two series overlap, permitting 
a cross-calibration of Wolf's index in terms of the 
aa-index. No magnetic data is available for 1805- 
1813 so Wolf used yearly counts of auroral displays 
and calibrated them in terms of magnetic activity. 
Finally, Wolf wanted to emphasize the sunspot cycle 
variation of geomagnetic activity and had smoothed 
his index to remove short-term fluctuations (though 
he ends up reporting monthly values!). 
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Since the year 1900 the general level of geomagnetic 
activity has increased at least two-fold even to the 
extent that the activity at recent sunspot minima ex- 
ceeds that of the sunspot maximum of 1906. It seems 
•likely that either the average solar wind speed and/ 
or the magnitude of the interplanetary magnetic 
field has changed significantly over the past 75 
years. A general change in solar wind speed entails 
a corresponding change in field magnitude because of 
different spiral angle of the interplanetary magnetic 
field. The product BV is in fact given' by 

BV = B R (V 2 +n 2 R 2 ) % (21) 

where B R is the radial component of the field at dis- 
tance R and ft = 2.87 yradian/s is the angular velocity 
of solar rotation. If V becomes considerably smaller 
than ftR = 430 km/s at 1 AU, the quantity BV is almost 
constant provided that B^. is constant. The net result 
is that an overall change in solar wind speed alone 
changes the am-index according to am'vV rather than 
am^V 2 , which means that the decrease by a factor of 
two in am-index back to 1900 corresponds to a similar 
decrease of V from 400 km/s to 200 km/s. It seems 
more likely to the present author that it is mainly 
the field magnitude that has changed through a change 
of B r «... a change of the open magnetic flux). 

The sunspot number is presumable some crude measure 
of the closed magnetic flux and it it not unreason- 
able to expect -a similar trend in the open flux. As 
seen in Figure 24 the long-term trends in geomagnetic 
activity since 1781 follows rather closely similar 
trends in the maximum sunspot numbers in each cycle, 
suggesting similar trends in available magnetic flux. 
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It would be very desirable to be able to infer the- 
interplanetary magnetic field strength independently 
from the solar wind speed from geomagnetic data. Some 
recent works (Russell and Fleming, 1976; Gul'yel'mi 
and Bol'shakova, 1973) suggest that the period of cer- 
tain types of micropulsations of the geomagnetic 
field (Pc2 to Pc4 with periods from 5 to 150 seconds) 
is strongly controlled by the magnitude of the inter- 
planetary magnetic field and can in fact be used as a 
diagnostic for that quantity. No attempt has yet 
been made to extend the analysis to the pre-satellite 
era . 


11. CONCLUSION 

In this chapter we have outlined some current ideas 
about the interaction between the solar wind and the 
earth's magnetosphere and showed how they can help or- 
ganize the analysis of how solar wind parameters in- 
fluence geomagnetic activity. A rather detailed de- 
scription of the dependence of geomagnetic activity 
on solar wind conditions is now available and permits 
conclusions about large-scale and long-term proper- 
ties of the sun and of the solar wind to be deduced 
from the geomagnetic records. On the other hand, the 
discovery of coronal holes seems to have brought us 

close to the solution of the old problem of a solar 

*\ 

identification of the elusive M-regions though'f re- 
sponsible for recurrent sequences of enhanced geomag- 
netic activity. Thus solar and geophysical research 
have recently strengthened their links in the impor- 
tant study of the environment of mankind. 
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